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ABSTRACT 

Observations of local star forming galaxies have revealed a correlation between the rate 
at which galaxies form stars and their X-Ray luminosity. We combine this correlation 
with the most recent observational constraints on the integrated star formation rate 
density, and find that star forming galaxies account for 5-20% of the total soft and 
hard X-ray backgrounds, where the precise number depends on the energy band and 
the assumed average X-ray spectral energy distribution of the galaxies below ~ 20 
keV. If we combine the Lx-SFR relation with recently derived star formation rate 
function, then we find that star forming galaxies whose X-ray flux falls well (more 
than a factor of 10) below the detection thresholds of the Chandra Deep Fields, can 
fully account for the unresolved soft X-ray background, which corresponds to ~ 6% 
of its total. Motivated by this result, we put limits on the allowed redshift evolution 
of the parameter cx = ^x/SFR, and/or its evolution towards lower and higher star 
formation rates. If we parametrize the redshift evolution of cx c>c (1 -t- z)^, then we find 
that 6 ^ 1.3 (95% CL). On the other hand, the observed X-ray luminosity functions 
(XLFs) of star forming galaxies indicate that cx may be increasing towards higher 
redshifts and/or higher star formation rates at levels that are consistent with the 
X-ray background, but possibly at odds with the locally observed ix-SFR relation. 

Key v^rords: galaxies: high redshift - galaxies: stellar content - X-rays: binaries - 
X-rays: galaxies 



1 INTRODUCTION 

Galaxies contain various sources of X-ray emission which in- 
clude: (i) active galactic nuclei (AGN), which are powered 
by accretion of gas onto a supermassive black hole (ii) hot 
(T ^10® K) interstellar gas, and {iii) X-ray binaries, which 
consist of a compact object, either a neutron star or a stellar 
mass black hole, and a companion star from which the com- 
pact object accretes mass. The X-Ray luminosity of star 
forming galaxies without AGN is dominated by so-called 
high-mass X-ray binaries (HMXBs, e.g. Grimm et al. 2003), 
in which the neutron star or stellar mass black hole is accret- 
ing gas from a companion that is more massive than ~ 5iV/o . 
HMXBs are thus tightly linked to massive stars, and since 
massive stars are short lived, the combined X-ray luminosity 
of HMXBs is expected to be Unked to the rate at which stars 
form (e.g. iHelfand fc Moranll200ll ). The X-Ray luminosity 
of star forming galaxies is indeed observed to be correlated 
with the rate at which they are forming stars (|Grimm et al.l 



E-mail:dijkstra@mpa-garcliing. mpg.de 



2002 



2004_ 
l2011b 



Ranalli et al 



Mineo et al.l 



20031: iGilfanov et"all 12004 jPersic et al 



2011a: iLehmer et al.l [20101 : iMineo et a" 
ISvmeonidis et aflboill ) 



This 'Lx-SFR relation' encodes a wealth of informa- 
tion on various astrophysical processes. These include the 
initial mass function (IMF) of stars, the fraction of mas- 
sive stars that form in binaries, the mass ratio distribution 
of binary stars, the distribution of their separations, the 
gas metallicity, and t h e common envelope efficiency (e.g. 
iBelczvnski et af I I2OO2I : iMineo et all l2011bl . and references 
therein). Despite its dependence on many astrophysical pro- 
cesses, the Lx-SFR r elation is observed to hold over ~ 4 
orders of magnitude (|Mineo et al.l [2011bl ). with a modest 
scatter (cr = 0.4 dex). 

Existing observations have only been able to probe the 
Lx-SFR relation in nearby galaxies. There are observational 
hints-as well as theoretical expectations- that cx is higher 
in low mass galaxies and/or low metallicity enviro nments 
(Dravl l2006l : iLinden et a"l]|2010l : iKaaret et aLll201ll '). which 
suggests that cx = Lx/SFR co uld have been substa ntially 
higher at higher redshifts (see iMirabel et al.l I2OIII . for a 
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summary). However, quantitative constraints on ratio cx at 
higher redshifts are virtually non-existent. The main reason 
is that the X-ray flux that reaches Earth from individual 
star forming galaxies typically falls well below the detection 
threshold of existing X-ray telescopes (see §[2|. 

The relation between SFR and Lx also plays an 
important role in determining the thermal history of 
inter g alactic medium at ver y high redshifts (e.g. lOhl 
I2001I : IVenkatesan et al] l200ll ). and strongly affects the 
21-c m signal from atomic hydrogen during the dark 
ages (iFurlanetto e t al.' 200 61: iPritchard fc Furlanettd [ioOTi : 
IPritchard fc Loebli2010. ; Alvarez et alj|2010l '). Current theo- 
retical models of this reheating process explore values for cx 
that span ~ 4 orders of magnitude. 

The goal of this paper is to investigate whether it is 
possible to put any constraints on the evolution of Cx with 
either redshift and/or towards high/low star formation rates 
than probed by existing observations of individual galaxies, 
using the cosmic X-Ray backgrounds (CXBs). 

Our paper is organized as follows. In § [S] we summarize 
recent observational constraints on the levels of the total X- 
ray background (XRB, both soft and hard), as well as the 
fraction of the XRB that has been resolved in discrete X- 
ray sources. In §[3] we compute the total contribution of star 
forming galaxies to the soft and hard X-Ray backgrounds, 
and show that this contribution can be substantial. In § [3] 
we show that star forming galaxies, too faint to be detected 
as individual X-ray sources, can fully account for the un- 
resolved portion of the XRB. In §[S]we put constraints on 
possible trends in the Lx-SFR relation using the unresolved 
soft XRB. There are several additional candidate sources 
which contribute to the unresolved CXBs (these include for 
example low mass X-ray binaries and weak AGN, see § [6]). 
To constrain the evolution of cx we will allow the entire un- 
resolved SXB to be produced by HMXBs (and the hot ISM) 
in X-ray faint star-forming galaxies. This results in conserva- 
tive upper limits on possible evolution in cx. Finally, we con- 
clude in § El The cosmological parameters used throughout 
our discussion are (O ^, Ha, fe) = (0.27,0.73,0.046,0.70) 
l|Komatsu et al.ll2009l '). 



2 THE COSMIC X-RAY BACKGROUNDS 
(CXB) 

2.1 The Observed Soft X-ray Background (SXB) 



(£;=l-2 keV) CXB (SXB) 



amounts 
J .-2 



The total soft 

Si-2 = 4.6 ± 0.3 X 10 

iHickox fc Markevitc B I2OO6I '). iHickox fc MarkevitchI ((20oi 



to 



find an unresolved SXB intensity of 1.04 ± 0.14 x 10"^^ 
erg s~^ cm~^ deg~^, after removing all point and extended 
sources detected in the Chandra Deep Fields (CDFs). The 
detection threshold in the Chandra Deep Field-North (CDF- 
North) c orresponds to Sth ~ 2.4 x 10~^^ erg s~^ cm^^ 



lAlexander ct al, 2003; IHickox fc Markevitch|[2007a ') 



Hickox fc MarkevitchI l|2007al 'l showed through a stack- 
ing analysis that ~ 70% of the unresolved component is 
accounted for by sources that are detected with the Hubble 
Space Telescope, but not individually as X-ray sources. Their 
stacking analysis shows that these X-ray undetected sources 
have an average X-ray flux that is (s) ~ 0.15 — 0.30sth ~ 



3.6-7.2x10"^* ergs"^ cm~ MHickox fc MarkevitchI (|2007bl ) 
find that the cumulative number of these X-ray unde- 
tected HST sources brighter than some of X-ray fiux s, 
is well-described by a power-law A*'(> s) oc s"^ , where 
/3 = l.llg'g. We can then estimate the minimum X-ray 
fiux, Sniin, that the HST-detected sources probed from 
(s) = and find that = 

1.2 - 3.2 x 10"^* erg s"^ cm"^ (for /3 = 1.1). Through- 
out, we use the remaining 30% of the unresolved SXB as 
a measure of the 'true unresolved SX B', which amounts to 
3.4±1 .4 X 10"" erg s"^ cm"^ deg"^ l|Hickox fc MarkevftchI 
l2007al ). down to a minimum fiux Smin. 

2.2 The Observed Hard X-ray Background (HXB) 



The total hard (£=2-8 keV) CXB (HXB) amounts to 
S2-8 = 1-7 ± 0.2 X 10~ " erg s~^ cm"^ deg" ^ (e.g . 
iHickox fc MarkevitchI 120061 ). iHickox fc MarkevitchI (120061 ) 
find an unresolved HXB intensity of 3.4 ± 1.7 x 10~^^ erg 
s~^ cm~^ deg~^. We will not use the unresolved HXB to put 
constraints on the redshift evolution of cx for two reasons: 
(i) the uncertainties on the unresolved HXB are larger than 
for the SXB, and {ii) as we will explain below (§ I3.2|l we 
are much more sensitive to the assumed X-Ray spectrum 
outside the Chandra bands, when we compare our models 
to the HXB. 



3 THE CONTRIBUTION OF STAR FORMING 
GALAXIES TO THE CXB 

3.1 The Model 

The total contribution ^tot (in erg s^^ cm^^ deg"^) of star 
forming galaxies to the SXB is given by (Appendix iB|) 



An c 



dz 



{l + zY£{z) 



p,(z)Cx{z,V). 



(1) 



Here ASTi ~ 3.0 x lO"* sr deg-^ £{z) = + + 

and p*[z) denotes the comoving star formation rate density 
at redshift z (in Mq yr~^ cMpc""^, where 'cMpc' stands 
for c o-moving Mpc). We adopt the star formation history 
from iHopkins fc BeacomI (|2006l ). u sing the pa r amet ric form 
p*(0) = {a^bz)h/{lJr{z/cf) from iCole etaH (|200li ). where 
a = 0.017, 6 = 0.13, c = 3.3, and d = 5.3. We note that the 
amplitud e of the function Pt(z) d epends on the assumed 
IMF fsee lHopk ms fc Beacom|[200^ , for a more detailed dis- 
cussion). The ado pted normaliza tion derives from classi- 
cal Sa lpeter IMF (Salpeter '1955*) between 0.1 to IOOMq 
l|Hopkina.2004 ). The same IMF w as assumed in the d eriva- 
tion of the SFR-Lx relation (see iMineo et al.]|2011bl ). and 
our calculations are therefore self-consistent. 

The term jC{z, F) = cxK{z, F) denotes the 'K-corrected' 
X-ray luminosity (in erg s"'^) per unit star formation rate 
in the observed energy range E1-E2. In this paper, Ei = 1.0 
keV and E2 = 2.0 keV for the soft band, ar id Ei = 2.0 
keV an d E2 — 8.0 keV for the hard band. iMineo et al.l 
lj2011bl ) recently determinecQ the value of cx = Lx/SFR, 



^ The literature contains values for cx that at face value ap- 
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Figure 1. The top panels show the fraction of the total soft {left) and hard (right) X-Ray backgrounds that can be attributed to star 
forming galaxies, as a function of the assumed photon index F, for power- law X-Ray spectral energy distributions (see text). We have 
drawn curves for a range of observed values for cx = ix/SFR (see text). Depending on the assumed value for F, star forming galaxies 
can account for ~ 5 — 15% of the total observed soft X-ray background, and up to ~ 20% of the hard X-ray background. For comparison, 
ISwartz et al. I l l2004) found that typically ULX spectra in the Chandra bands could be described by a powerlaw with (F) = 1.7, which 
would place the contribution of star forming galaxies to the soft X-ray background at ~ 9 — 13%. As the spectra of ULXs at £" > 10 
keV are poorly known, the contribution of star forming galaxies to the hard X-ray background is more uncertain. The lower panels show 
that the contribution per unit redshift, dStot/d.z, peaks at low redshift z ~ 0.4 — 0.5. 



where SFR denotes the star formation rate in Mq yr~^, to 
be Cx = 2.6 X 10^^ erg /[Mq yr~^] when onl y compact re- 
solved X-Ray sources in galaxies are included. iMineo et al.l 
l|2011br ) also found that the best fit Cx.max = 3.7 x 10^^ 
erg s~^ [Mq yr~^]~^ for unresolved galaxies in the Chandra 
Deep Field North and ULIRGs. However, a non-negligible 
fraction of this additional unresolved flux is in a soft com- 
ponent, and would not contribute to soft X-Ray background 
(measured in the 1-2 keV band, see Bogdan & Gilfanov 
2011). We will not attempt to model in detail the contribu- 
tion of unresolved X-Ray emission. Instead, Figures[Tl[2l and 



pear both lower (e .g. [Persic fc RephaelH l2007tl and higher (e.g. 
iRanaUi et al]|2003f) by a factor of a few. However, some studies 
measured the X-Ray luminosity (Lx) in the range 2-10 keV (e.g . 
Gilfanov et"al] 12004 iPersic et al.ll2004 iPersic fc Rephaelill20oi 
Lehmer et al.ll2010l) , or 0.5-2.0 keV llRanalli et al.ll2003h . Further- 
more, some studies have derived values for cx using a forma- 
tion rate of stars (SFR ) in the mass range 5 M/Mq ^ 100 
(e.g. iPe rsic et al]|2004l 'l. which results in larger values for cx. 
iMineo et al.l l l2011bl) ^ discuss that most studies are consistent with 
their derived value when identical definitions for 'SFR' and '-Lx' 
are used. 



[4] show our results for the full range Cx = 2.6 — 3.7 x 10^^ erg 
s"^ [Mq yr"^]"\ We point out that iMineo et all (|2011bi ') 
measured the X-ray luminosity over the energy range 0.5 — 8 
keV. 

This observed relation between Lx and SFR holds over 
~ 4 orders of magnitude in SFR from SFR~ 0.1 - 1000 Mq 
yr^^, with a scatter of cr ~ 0.4 dex. We include this scatter 
in our calculation by convolving equation ([1} with a lognor- 
mal probability distribution function for cx with a mean of 
(l ogcx) = 39.4 and st andard deviation of o" = 0.4 (see Fig. 9 
of|M meo et al.ll2011bl ). Because the log-normal distribution 
is symmetric in logcx, it is skewed towards larger values of 
Cx in linear coordinates. The observed scatter in the Lx-SFR 
relation therefore enhances our computed contributions to 
the X-ray backgrounds. In the absence of this scatter, Stot 
reduces by a factor of exp (— ia^ln^lO) ~ 0.65. 

Finally, to compute the contribution of a galaxy at 
redshift z to the CXB in the observed energy range Ei- 
E2, we need to compute the galaxy's luminosity in the 
range [Ei-E2]x(l -I - z) keV. In ana logy to the standard 
'K-correction' (e.g. iHogg et al.l |2003 ). we multiply cx by 
Kx{z,r) = iIei{1 + z),E2(l + z))/l{0.5keV,8.0keY). 
Here, I{x, y) = JJ' En{E)dE, where n{E)dE denotes the 
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number of emitted of photons in the energy range E±dE/2. 
Note that we assume that the X-Ray SED does not depend 
on the star formation rate tp. In this paper, we explore power 
law SEDs for the form n{E) oc , and the integrals I{x, y) 
can be evaluated analytically. Finally, we take the redshift 
integral from Zmin = to Zmax = 10. The results are only 
weakly dependent on the integration limits we pick, as long 
as 2max ^ 4.0 (see the lower right panel of Fig. [2|. 



3.2 Results 

The top panels of Figure [T] show the fraction of the total 
soft {left panel, 1-2 keV) and hard [right panel, 2-8 keV) 
X-Ray backgrounds that can be attributed to star forming 
galaxies, as a function of the assumed photon index F, and 
for a realistic range of cx (see above). Depending on the 
assumed value for F, star forming galaxies can account for ~ 
5— 15% of the total soft X-ray background, and up to ~ 20% 
of the observed ha rd X-ray background. For comparison, 
ISwartz et all (|2004l ) found that the observed distribution of 
F for ultra luminous X-ray sources (ULXs) had a mean (F) = 
1.7, and mode (i.e. most likely value) of Fp^ ~ 2.0 in the 
Chandra bands. For the range F = 1.7—2.0 the contribution 
of star forming galaxies to the soft X-ray background is ~ 
8 — 13%. The contribution to the hard XRB depends more 
strongly on F because for this calculation, the K-correction 
involves a larger extrapolation of the assumed SED, which 
is poorly known at _E > 10 keV. 

The lower panels show the differential contribution as 
a function of redshift. This differential contribution is given 
by dStot/dz. These plots show that dStot/dz peaks at 2 ~ 
0.38, and that the dominant contribution to the total X- 
Ray background comes from lower redshifts. Indeed, ~ 50% 
[~ 25%] of the total contribution comes from 2 ^1.3 [z <;0.8] 
(see §112]). ^ 



4 CONTRIBUTION OF X-RAY FAINT STAR 
FORMING GALAXIES TO THE SXB 

4.1 The Model 

We can compute the contribution Sx of star forming galax- 
ies, fainter than some observed soft X-ray flux Smax, to the 
SXB as 



Sx 



AO, 



dz 



f 

Jo 



47tHoJo il + z)^£(z) 
dlogLx n(logLx,2) LxKx{z,T), 



(2) 



where n(log Lx, z)dlog Lx denotes the comoving number 
density of star forming galaxies with X-ray luminosities in 
the range logI/x±oilogLx/2 (i.e. the units of n(logLx, z) are 
cMpc~^ dex~^). Here, Lx denotes the X-ray luminosity of 
galaxies in the 0.5-8.0 keV (restframe). The integral over Lx 
then extends up to Lx.max = 4:nd1{z)snia.x/Kx{z,T), where 
di,{z) is the luminosity distance to redshift 2. The quantity 
n(logLx), also referred to as the X-ray luminosity function 
(XLF) of star forming galaxies, is given by 



where n{%l),z)d^l) denotes 'star formation rate function', 
which gives the comoving number density of galaxies that 
are forming stars at a rate SFR= il)±dtp/2 at redshift 2. The 
function P(log Lx|V')'^log Lx denotes the probability that a 
galaxy that is forming stars at a rate ip has an X-ray lu- 
minosity in the range logLx ± dlogLx/2. We describe both 
functions in more details below. We start the integral over 
at 7/)miii = 1O~^M0 yr~^, which corresponds approximately 
to the SFR that is theoretically expecte d to occur in dark 
matter halos of mass Mhaio ~ IO^M r l|Wise fc Cenll2009l : 
iTrac fc Ccn '2007l: IZheng et al.llioiol 'l. Our final results de- 
pend only weakly on V'min (see the upper left panel of Fig. (2]). 
The ?/;-integral extends up to i/imax ~ 10^ yr~^, with our 
results being insensitive to this choice. 

In the local Universe (2 ~ 0), the function n{tl;, z) [units 
are cMpc"'^ (A^© yr~^)~^l ap pear to be describe d accu- 
rately by a Schechter function (jBothwell et al.|[201ll ) 

1^(|^)"'-'--' w 

where a = -1.51 ± 0.08, $* = (1.6 ± 0.4) x 10"" cMpc"^ 
and t/j* = 9.2 ± O.SM© yr~\ The redshift evolution of 
n(i/), 2) is not well known. We assume throughout that 
a{z) = —1.51 — 0.23(?(2), where g{x) = -arctan a; is a func- 
tion that obeys g{Q) = and lim g{z) — 1. This steepening 

2—^00 

of the low-end of the star formation rate function at higher 
redshifts reflects the steepening of UV luminosity functions 
towards higher redshifts (e.g. Arnouts et al. 2005, Reddy 
& Steidel 2009, Bouwens et al. 2006, 2007,2008), and the 
observation that dust-obscuration is negligible for the UV- 
faint galaxies (e.g. Bouwens et al. 2009). The factor '-0.23' 
causes a — ^ —1.74 at high redshift, which corresponds to 
the best-fit slope of the UV-luminosity function of 2 = 6 
drop-out galaxies (Bouwens et al. 2007). The redshift evolu- 
tion of $* and -i/)* is more difficult to infer from the redshift 
evolution of the UV luminosity functions, because of dust. 
We have taken two approaches: we constrain either the red- 
shift evolution of $* or T/i'-while keeping the other fixed- to 
match the inferred redshift evolution of the star formation 
rate density (see Appendix [X] for more details). In reality, 
we expect both parameters to evolve with redshift, and that 
our two models bracket the range of plausible more realistic 
models. 

The function P{\og Lx\i^) is given by a lognormal dis- 
tribution 



P(logLx!^) = 



2na 



exp 



log 



(-c-x>> 



2cr2 



(5) 



n(logLx, z) = 



d^ n{^p, z)P {log Lxlip), 



(3) 



where (Lx) — cx x denotes the X-ray luminosity (mea- 
sured in the 0.5-8.0 keV rest frame) that is expected from the 
observed SFR-Lx relation. The standard dev iation a = 0.4 
deno t es the observed scat t er in this relat i on (|L"ehmer et al.l 
I2OIOI : iMineo et al.|[2011bl ). iRanalli et all l|2005l ') made very 
similar predictions for the XLFs, but instead of using 
star formation rate functions, t hey used gala xy luminos- 
ity functions in various bands. iRanalli et al.l 1^2005,) also 

^ Previous work showed that this star for mation rate functio n 
can be described by a log-normal function l lMartin et al.|[2005^ . 
However, this lognormal function does not provide a good fit to 
the observed star formation rate function particularly for low and 
high star formation rates (see Fig 4 of Bothwell et al. 2011). 
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Figure 2. The contribution to the soft X-Ray background (SXB, E=l-2 keV in the observer's frame) by galaxies whose individual 
soft X-ray flux is less than Smax- The median unresolved SXB is represe nted by the black solid horiz ontal lines, and its 68% confidence 
levels by the the gray region, bounded by black dashed lines (taken from lHickox &: Markevitchll2007al ).. The red solid lines show Sx as a 



function of various model parameters. Our fiducial model assumes F = 2.0, Sn 



■ 2.4 X 10-1* ergs-l 



: O.OIM, 



, and 



.^max = 10. The upper left panel shows that i/imin only weakly affects Sx, because the faint end of the star formation function (especially 
at low z) is not steep. The upper right panel shows that Sx depends weakly on F, provided that F <;2. The lower left panel shows that Sx 
also depends weakly on Smax, unless Smax ^10—^* erg s" ^ cm"^. The dotted vertical line shows the X-ray detection threshold in CDF-N 
l|Alexander et al.ll20o3 ). The black filled circle on this line shows the unresolved SXB derived by Hickox & Markevitch (2006, i.e. before 
subtracting the contribution from X-Ray faint HST detected sources). The gray region here brackets the effective minimum X-ray flux 
Sjnin that is probed by stacking X-ray undetected HST sources (see ^ 12.11 1. The lower right panel shows that Sx again depends weakly 
on 2max, provided that Zmax > 4. See the main text for a more detailed interpretation of these plots. These plots show that star forming 
galaxies that are too faint to be detected as individual X-ray sources, can account for the full unresolved SXB, an d that this stateme nt 
is insensitive to details in the model when F <2, which is reasonable given the available observational constraints jSwartz et aI.ll2004^ . 



assumed a lognormal conditional probability functions for 
P(logI/x|iY), where Ly denotes the galaxy luminosity in 
some other band Y. 



4.2 Results 

Our fiducial model assumes F — 2.0 (which corresponds to 
Fpk, see above), Smax ~ 2.4 x 10"^* erg s-^ cm"^ (close to 
the middle of the range for Smin that was quoted in § 12. 1|) 
, 'i/'min = W~^Mq yr~i, and Zmax ~ 10. As mentioned pre- 
viously, we explore two choices for extrapolating the star 
formation rate function with redshift, which likely bracket 
the range of physically plausible models. When we evolve 
tp*{z), but keep $* fixed, we find that our fiducial model 
gives Sx = 2.4 x 10"^^ erg s"^ deg~^ cm"^. On the other 
hand, when we evolve ^*{z), but keep t/i* fixed, we obtain 
Sx = 2.6 X 10~" erg s"^ deg"^ cm"^. The fact that this 
difference is small is encouraging, and suggests that our ig- 
norance of the star formation rate function at z > does 



not introduce major uncertainties into our calculations. We 
have verified that the second model generally yields slightly 
higher values for Sx- To be conservative, we focus on the 
first model in the reminder of this paper. 

Figure [2] has four panels, each of which shows the me- 
dian unresolved soft X-ray background {black solid lines), 
and its 68% confidence levels (the gray region, bounded by 
black dashed lines) from Hickox & Markevich (2007b). The 
red solid bands show the contribution from galaxies whose 
individual soft X-ray flux (1-2 keV observed frame) is less 

than Smax- 

• In the upper left panel we plot Sx as a function of 
V'min- We find that 5x depends only weakly on ipmin- That 
is, very faint galaxies do not contribute significantly to 
Sx- This is because the majority of the contribution to Sx 
comes from galaxies at z < 2 (see below, and Fig. [TJ , where 
the 'faint' end of the star formation function increases as 
oc tli~'^'^, and the overall star formation rate density is domi- 



© 2009 RAS, MNRAS OOO.IlM 



6 Dijkstra et al. 



nated by galaxies that are forming stars at a rate close to -i/)*. 

• The upper right panel shows 5x as a function of F. 
We find that Sx decreases with F. As most of the contri- 
bution to S'x comes from galaxies at z < 2, we are most 
sensitive to the X-ray emissivity of star forming galaxies at 
S = [1 - 2] X (H- 2) < [3 - 6] keV (restframe). For fixed 
Lx, increasing F reduces the fraction of the emitted flux 
at these 'higher' energies for steeper spectra (i.e. most of 
the energy lies near E — 0.5 keV for the steepest spectra), 
which results in a decrease in Sx- 

• The lower left panel shows Sx as a function of Smax. 
The vertical dotted line shows t he detection threshold 
in the Chandra Deep Field-North (|Alexander et al.1 12OO j : 
iHickox fc Markevitchlr2007ah . The gray region here brackets 
the effective minimum X-ray flux Smin that is probed by 
stacking X-ray undetected HST sources (see §[2TT]). Clearly, 
Sx depends only weakly on Smax, unless Smax ;^10~^* 
erg s~^ cm~^. This weak dependence on Smax at larger 
fluxes can be easily understood: most of the contribution 
to Sx comes from z <2 (see below). For galaxies at 

corresponds to 
^, which 



-18 



cm 

ri40 



z = 1, Smax = 2.4 X 10 " erg s 
Lx = 47rd£(2)smax/i:x(2,F) = 2.9 X 10*" erg s , 
requires = 9.7Mq yr~^, which is close to 1/)*. We are 
therefore practically integrating over the full UV-luminosity 
function. Boosting Smax therefore barely increases Sx 
further. 

• The lower right panel shows Sx as a function of Znmx. 
This plots shows that Sx evolves most up to Zmax ~ 2, and 
barely when Zmax >4. That is, galaxies at higher redshift 
barely contribute to Sx as we also showed in Figure[T] (unless 
the conversion factor cx between Lx and SFR changes with 
redshift, see § 

The main point of Figure [2] is that it shows that star 
forming galaxies that are too faint to be detected as individ- 
ual X-ray sources, can account for the full unresolved SXB, 
and that this statement is insensitive to details in the model 
when F <2 over the energy range 1-6 keV. This last re- 
quirement is very reasonable, given t hat th e mean observed 
(F) = 1.7 for the majority of ULXs (|Swartz et al.,,2004 '). 



5 CONSTRAINING THE REDSHIFT 

EVOLUTION OF THE SFR-Lx RELATION 

5.1 Constraints from the SXB 

Since our fiducial model already saturates the unresolved 
SXB, we can ask what constraints we can set on the redshift- 
dependence of the Lx-SFR relatior[f|. Clearly, the SXB can 
only put constraints on models in which cx increases with 



As we already stated in §[T] there are several additional candi- 
date sources which contribute to the unresolved SXB (see §[6]for a 
summary) . To constrain the evolution of cx we will conservatively 
allow the entire unresolved SXB to be produced by HMXBs (and 
the hot ISM) in X-ray faint star-forming galaxies. Obviously, this 
results in upper limits on possible evolution in cx- 




2x10" 
A (erg/s/[M 



5x10" 

/yr]) 



10 



40 



Figure 3. The unresolved SXB constraints on the parameters A 
and b for a redshift evolution parametrization of the form cx = 
= j4(1 -|- z)'. Models that lie in the red region saturate the 
unresolved SXB at > 3cr (see §[5TTJ- T he grey regio n denot es the 
value for cxiz = 0) = A derived by iMineo et al.l l|2011lj ) (the 
solid vertical line denotes their best-fit value). For this value of 
A, b <1.3 (at 2(t). 



redshift. We consider models for which 
Lx 



SFR 



cx = A{l + z)'' 



(6) 



with 6^0, and investigate the constraints that the SXB 
places on the parameters A and b. We compute Sx (see 
Eq. |3]) on a grid of models which cover a range of A and 
b. Figure [3] shows how many a{= 1.4 x 10~^^ erg s~^ 
cm^'^ deg~^) above the median observed unresolved SXB, 
Sobs = 3.4 X 10~^^ erg s"^ cm~^ deg~^, the models lie. 
Models that lie above the uppermost solid line, indicated 
by the label '3-ct' result in Sx — Sobs > 3(7. That is, these 
models significantly overproduce the unresolved SXB and 
are practically ruled out. The light grey region bounded by 
the dotted vertical lines indicat es the 68% confidenc e region 
for the observed value for cx bv lMineo et al.l l|2011bl ). where 
we assumed 0.4 dex uncertainty on cx. Their best-fit value 
for Cx is indicated by the solid vertical line. Figure [3] shows 
that for this value of A, b ^1.6 at ^Sa. If we marginaliz^B 
over A, then we also get b ^1.6. 

We also investigate changes in the SFR-Lx relation of 
the form 

( model I SFR < x Mq yr"^; 

Cx -> ,f X Cx < model II SFR ^ y Mq yr'^ (7) 

[ models III&IV z > Zc 

The upper left panel of Figure 2] shows Sx as a function of / 
for model I for x = 0.1 {red band) and x — 0.5 {green band). 
This plot shows that / ^3 — 5 for x = 0.1 is at odds with the 



■* We obtain this 
^ dA U{b\A)Pp,i, 



It 



marginalized upper limit U{b) = 
i{A). Here, }J{b\A) denotes the upper 
and P„ 



limit on b given A, and Pprior(^) denotes our prior on the 
probability density function for A. Since A = cx{z = 0), 
we took Pprior(^) to be a logn ormal distribution with 
(logcx) = 39.4, and a = 0.4 jMineo et al.l l2011bl) . Fi- 



nally, 

^max 



we adopted Amin = 10^*'^ 
= 10''0-'5 erg s-1[Mq yr'^J-^ 



erg s 



-l[M0 yr- 



and 
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Figure 4. Same as Figure [2] but we modify our fiducial model such that the X-ray emission from star forming galaxies is boosted by a 
factor of cx — > / X cx (i.e. Lx — > / x ^x) when SFR< x Mq yr~^ {top left panel), SFR> y Mq yr^^ (top right panel), and z> Zc {lower 
left panel) as a function of /. The lower right panel shows a model where we boost cx by a factor of 100 {upper blue band), 10 {middle 
red band) and 3 {lower green band) at redshift Zj^^p as a function of zjump- This figure shows that even small boosts (i.e. / ~ a few) 
for SFR< O.IMq yr~^ or 2 ^ 2 violates the constraints posed by the unresolved SXB. It also shows that boosting cx at the large SFR 
end barely affects Sx- While large 'jumps' in cx are not allowed if these occur at low redshift (i.e. z ^3), the unresolved SXB cannot 
constrain large jumps at z > 5. 



unresolved SXB. This can be understood from the top left 
panel of Figure [21 which shows that adopting -i/imin = O.IMq 
yr"'^ results in Sx ~ 1-9 x 10"^"^ erg s~^ cm~^ deg~^, which 
corresponds to a reduction of ~ 25%. Hence, galaxies with 
< 0.1 contribute ~ 25% to S'x (for / = 1). 



0.001< „''^ _i 

Mq yr 

Boosting their contribution by a factor / >3 — 5 (the exact 
number depends on the precise fiducial choice for cx) causes 
Sx to exceed the unresolved SXB. 

The upper right panel shows that the soft XRB only 
allows constraints to be put on / ^10, and only if y <,25. 
Strong constraints on / are not possible for large ip. This is 
because galaxies that are forming stars at a rate ip ^IOOMq 
yr~^ are deep in the exponential tail of the star formation 
function. As a result of their small number density, they 
barely contribute anything to Sx- 

The lower left panel shows that boosting cx at 2 ^ 2 
by factors greater than / ^3 — 5 is again at odds with the 
unresolved SXB. This is because galaxies at a > 2 contribute 
noticeably to 5x for our fiducial choice of cx. However, the 
unresolved SXB cannot place tight limits (yet) on cx at very 
high redshifts. This is illustrated in the lower right panel 
where we show S'x for model IV: once zjump ^5, boosting cx 
by as much as ~ 100 has little impact on Sx. 



5.2 Constraints from the Galaxy XLFs? 

As part of our analysis, we compute theoretical galaxy 
X-Ray luminosity functions (XLFs) using equation 
Ob served galaxy XLFs h a ve been presented, for e x ample , 



bv iNorman et al 



iGeorgakakis et al 



(200S 



(2006) 



(2004): iGeorgantopoulos et all ll2005h: 



Tzanavaris & GcorEantopoulo^ 
Tzanavaris fc Georga ntopoulos (2008) present 
galaxy XLF for late-type (i.e. star forming) galaxies in 
two redshifts bins: the first bin contains galaxies with 
< 2; < 0.4 (zmcd = 0.14), and the second bin contains 
galaxies with 0.4 ^ z < 1.4 (zmed = 0.68). Note that 
ITzanavaris fc Georgantopoulos! (|2008l ) derive X-Ray lumi- 
nosities in the 0.5-2.0 keV band, and we properly K-correct 
our models into their band. 

In Figure [5] the data points show the observed X-Ray 
luminosity functions (XLFs) from Tzanavaris fc Georgan- 
topoulos (2008), while the red dotted line shows the best- 
fit Schechter function derived by Georgakakis et al (2006). 
The solid lines show predictions for our fiducial model 
Cx = A{l + zf withyl = 2.6x10^^^ ergs"^ [Mq yr~^]~'^ , and 
b = 0.0. At z 0.4 our model provides a good fit at the two 
brightest X-ray luminosities, but overpredicts the number 
density of fainter sources by a factor of ~ 2 — 3. We obtain 
a better fit if we lower A = 1.4 x 10^^ erg s'^ /[MQyr~\ 
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which corresponds to ~ 0.2 dex , and thu s hes within the dis- 
persion that was found by Mineo et al.l jioilb ). This lower 
value of A corresponds to almost exactly the value quoted 
by Lehmer et al (2010, their j3, although these authors mea- 
sured Lx in the 2.0-10.0 keV range). At higher redshift, 
our model significantly underpredicts the number density at 
ix, 0.5-2.0 ^10*^ erg s^^. This may suggest that either that 
the XLFs (strongly) favor cx to increase towards higher red- 
shift, and/or SFR ^ a few tens of Mq yr~^ (see below). 
Alternatively, the observed XLFs of star forming galaxies are 
contaminated by low luminosity AGN, which are difficult to 
identify at these X-ray luminosities. 



Our p redicted XLFs agree quite well with previous pre- 
dictions bv lRanalli et al.l l|2005l ). for the model in which they 
assume that the redshift evolution of the luminosity func- 
tions is solely the result of evolution in the number density 
of galaxies (this is referred to as 'density evolution'). Their 
model also underpredicts the numb er density of lu minous 
X-ray sources at higher redshift. Ra nalli et al.l (|2005 ) found 
that a better fit to the high-redshift data is obtained for 
a model in which solely the luminosity of galaxies evolves 
('luminosity evolution') as oc (1 + zf-'^-'^-^. Our work also 
indicates that evolution in the number density of star form- 
ing galaxies is not enough to explain the observed redshift 
evolution of XLFs, and that some additional luminosity evo- 
lution is preferred. 



It is possible to compute the likelihood £(^4, &) = 
exp[— 0.5x^1 by fitting to the observed XLFs for any com- 
bination of A and 6 describing the redshift evolution of cx 
(see Eq. However, we found that this formal fit is dom- 
inated by the two lowest luminosity data points at z ^ 0.4, 
which lie significan tly be low the XLF that was derived by 
iGeor gakakis ct al.^ (|2006l ). Furthermore, we found that this 
procedure also depends somewhat on the assumed redshift 
evolution of the star formation rate function. We therefore 
do not pursue a detailed statistical analysis on constrain- 
ing the redshift evolution of cx with the XLFs. Instead, 
we show in Figure [6] an example of a model where cx in- 
creases both with redshift and at high SFR: this model with 
^ = 1.4 X 10^" erg s"^ [Mq yx'^Y^ , b = 1.0, and cx -> 3cx 
at SFR;^ 15M0 yr~^ fits the observed XLFs much better. 
The value 6 = 1.0 is consistent with the X-ray background 
constraint (see Fig. [3]), and boosting cx by a factor of 3 
at SFR^ 15Mq yr~^ is also consistent with the soft XRB 
(see the top right panel of Fig. 2)). However, boosting cx 
by a factor of 3 at SFR^ 15Mq yr~^ appears inconsistent 
with direct constraints on the Lx-SFR relation. In Figure [6] 
we used cx = 4.2(1 + z) x 10^^ erg s"^ [Mq yr"^]"^ at 
SFR ^15A/o yr"^ which translates to cx = 8.4 x 10^^ erg 
s"^ [Mq yr"^]~^ at z = l, while iMineo et al.l (|2011b l) found 
Cx = 3.7 X lO^'' erg s"^ [Mq yr"^]"^ for their sample of 
unresolved, high SFR, sources at 2; ~ 0.2 — 1.2. However, 
the sample o f high-S FR, high-z galaxies that was studied by 
[m ineo et aU (|2011bl ) is rather limited, and if the dispersion 
around this mean quantity is also 0.4 dex, then our model 
may be consistent with the observed dispersion around this 
value out to z ~ 1. 




Figure 5. The observed X-Ray luminosity functions (XLFs) of 
galaxies from Tzanavaris &: Georgantopoulos (2008, data points) 
and Georgakakis et al (2006, the red dotted line shows their best- 
fit Schechter function) are compared to our fiducial model cx = 
A{l + z)'' with A = 2.6 X 10^^ erg s"! [Mq yr"!]"!, and b = 0.0. 
At 2 ^ 0.4 our model provides a good fit at the two brightest X- 
ray luminosities, but overpredicts the number density of fainter 
sources by a factor of ~ 2 — 3. The z ^ 0.4 XLF prefers Cx to be 
lower b y ~ 0.2 dex, within the dispersion found by Mi neo et al.l 
J2011bl) . At higher redshift, our model significantly underpredicts 
the number density at Lx, 0.5— 2.0 ^10*^ ^rg s~^. This suggests 
that either fe > 0, or that cx increases at high ip (see text). 




log \.0.5-2 

Figure 6. Same as Figure [S] except that here we improved the 
fit to the observations by forcing cx to evolve with redshift and 
star formation rate. The model now assumes cx = A{l + z)'' with 
A = 1.4 X 10^^ erg s"^ [Mq yr~l]~l, and b = 1.0. Furthermore, 
we boosted cx by a factor of ~ 3 at SFR ^15Mq yr~^. This plot 
illustrates that the XLFs favor an increase of cx with redshift 
and/or at larger SFR. 



6 DISCUSSION & CONCLUSIONS 

Observations have established that a correlation exists be- 
tween the star formation rate of galaxies and their X-ray lu- 
minosity (measured over the range E=0.5-8 keV, e.g. Ranalli 
et al. 2003, Grimm et al. 2003, Gilfanov et al. 2004, Lehmer 
et al. 2010, Mineo et al. 2010). This 'Lx-SFR relation' en- 
codes a wealth of information on various astrophysical pro- 
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cesses, and strongly affects the thermal evolution of the in- 
tergalactic medium during the early stages of the epoch of 
reionization. Existing observations have only been able to 
probe this relation in nearby galaxies, and while theoreti- 
cally there are good reasons to suspect that Cx = Lx/SFR 
increases towards higher redshifts, observational constraints 
are virtually non-existent. 

In this paper, we have investigated whether it is possi- 
ble to put any constraints on the evolution of Cx with either 
redshift and/or towards high/low star formation rates than 
probed by existing observations of individual galaxies. As 
part of our analysis, we have computed that the observed 
'local' relation, when combined with the most observational 
constraints on the redshift-evolution of the star formation 
rate density of our Universe, implies that star forming galax- 
ies contribute ~ 5 — 15% of both the soft and ~ 1 — 20% of 
the hard X-ray backgrounds (see § [Sjl . The ranges that we 
quoted are for a range of photon index 1 < F < 3. The ob- 
served r of ULX spectra in the Chandra bands is described 
by a distri bution with a mean of (F) = 1.7, and a mode of 
Fpk ~ 2.0 (|Swartz et al.1 120041 '). For the range F = 1.7 - 2.0 
the fractional contribution of star forming galaxies to the 
soft X-ray background is ~ 8 — 13%. The contribution to the 
HXB remains uncertain, mostly because of a more uncertain 
K-correction at the correspond ing high energies (also see 
iTrever fc Lahavl [l996.; .Nataraian fc Almainil |2000| . for ear- 
lier calculations of the contribution of star forming galaxies 
to the CXBs). 

We have then taken the most recent observational con- 
straints on the star formation rate function, which gives the 
comoving number density of star forming galaxies as a func- 
tion of their star formation rate [denoted by n{ip,z)], and 
computed what the contribution of 'X-ray faint' star forming 
galaxies to the soft X-ray background (SXB, corresponding 
to 1-2 keV in the observed frame) is. We found that galaxies 
whose individual observed flux is s ^ Smax ~ 2.4 x lO"'^* erg 
s~^ cm~^ between 1-2 keV, i.e. more than an order of mag- 
nitude fainter than the detection threshold in the Chandra 
Deep Field-North (see ^ 12. ip . can fully account for the un- 
resolved SXB. This statement is insensitive to details in the 
model as longs as the photon index, averaged over the entire 
population of X-ray emitting star forming galaxies, is F <.2, 
which corresponds to a very reasonable range given the ex- 
isting observational constraints on this parameter (§|31). 

Motivated by our result that X-ray faint star form- 
ing galaxies can fully account for the unresolved SXB, we 
put constraints on the redshift evolution of the param- 
eter Cx. When we parametrize the redshift evolution as 
Cx ~ A{1 + zf" , we found that the unresolved SXB requires 
that 6 ;^1.6 (3o"). We have also ruled out models in which cx 
is boosted by a factor of / ^2 — 5atzjsl — 2 and/or SFR 
^ 0.1 — O.5M0 yr~^, as they overproduce the unresolved 
SXB (see left panels of Fig. [J). We have found indications 
in the observed X-ray luminosity functions (XLFs) of star 
forming galaxies that cx is increasing towards higher red- 
shifts and /or higher star formation rates , but caution that 
this may indicate the presence of unidentified low luminosity 
AGN. The unresolved SXB allows for larger changes in cx 
at large values for SFR (see the top right panel of Fig. |4]). 
Finally, we also found that the SXB puts weak constraints 



on possible strong evolution (/ ~ 100) at 2 > 5 (see the 
lower right panel of Fig. [4](f|. 

There are many other undetected candidate sources 
which likely also contribute to th e unresolved SXB. These 
are briefly summarized below (see lDiikstra et al]|2004l . for a 
more detailed summary): 

• Observed AGN account for ~ 80% of the SXB. It would 
be highly unlikely that fainter AGN- i.e. those are too faint 
to be detected as individual X-ray sources-do not provide a 
signiflcant contribution to the unresolved SXB. 

• Our attention has focused on HMXBs, but low mass 
X-Ray binaries (LMXBs)-in which the primary has a mass 
<;5M0, dominate the X-ray luminosity of galaxies for 
which the speciflc star formation is sS FR <;10~^'' yi'~^ 
(|Gilfanov et al.ll2004l : ILehmer et aLlbOlO j). LMXBs give rise 
to a correlation between X-ray luminosity and tot al stellar 
mass , M,, which is Lx.lm xb ~ 9 x 10^* M, erg s~^ l|Gilfanovl 
l2004l : lLehmer et al.|[2010l '). In Appendix[C]we repeat the cal- 
culation of § |4] and replace the star formation rate function 
n(i/), z) with the stellar mass function n(A/*, z), and appro- 
priately replace P(logLxlV') with P(log LxlAf*). We found 
that faint 'quiescent' galaxies contribute about an order of 
magnitude less to the SXB than faint star-forming galaxies. 

• Thomson scattering of X-rays emitted predominantly 
by high-redshift sources can cause 1.0-1.7% of the SXB to be 
in a truly diffuse form (Soltan 2003 ). Similarly, intergalactic 
dust could scatter X-rays by small angles into diffuse halos 
that are t oo faint to be det ected individually (|Petric et al.l 
2006; Dii kstra fc L"oebll2009l ). 

• Wu fc Xud (2001) computed that clusters and groups of 
galaxies possibly contribute as much as ~ 10% of the total 
SXB. 

• A (hypothetical) population of 'miniquasars' powered 
by intermediate mass black holes m ay have contr i buted 
to ionizing and heatin g the IGM (|Madau et al.l |2004| : 
iRicotti fc OstrikeJ 120041 '). These miniquasars would emit 
hard X-ray photons that could cont ribute signiflcantly to the 
soft and hard X-ray backgrounds (|Ricotti fc Ostrikeill2004l : 
iDiikstra et"all l2004') . 

The likely existence of these additional contributors to 
the unresolved SXB implies that our constraints are conser- 
vative, and that actual limits on the redshift evolution of cx 
should be tighter. 

After our paper was submitted, Cowie et al. (2011) com- 
pared the average X-ray fluxes (obtained by a stacking anal- 
ysis) and restframe UV flux densities of sources with known 
redshifts in the 4 Ms exposure of the CDF-S field. Cowie 
et al. (2011) showed that this ratio -after an extinction 



5 iTreister et"al] 1I2OIII ') stacked 197 HST detected candidate 2 ~ 6 
galaxies and found significant X-ray detections in both the soft 
(0.5-2.0 keV) and hard (2.0-8.0 keV) bands (but sec Cowie et 
al. 2011 and Willott 2011). They derive an average rest frame 
2-10 keV luminosity of Lx,2-io = 6.8 X lO*'^ erg s~^, which they 
associate with obscured AGN. We can use this detection to place 
an upper limit on cx at 2 = 6. The mean star formation rate 
-averaged over the UV-luminosity function in the range —21.5 < 
A'/uv < —18.0, and not corrected for dust— is ~ 2Mq yr~^. 
The stacked X-ray detection therefore puts an upper limit on the 
boost / ^'X,2-io/l-5cx ~ 10'^. We verified that such a boost 
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correction- was consistent with the local Lx-SFR relation 
up to 2 ~ 4. The stacking of many source allowed Cowie et 
al. (2011) to probe down to s ~ Smax/4 = 5 x 10~^^ erg s"^ 
cm~^, which transla tes to a luminosity o f Lx ~ 0.3 — 4 x 10^" 
erg s~^ at z — 1 — 3. ICowie et al.l (|201lh therefore probe the 
redshift-evolution of cx at SFR ^1 — IOM0 yr~^ at these 
redshifts (depending on z and cx). For comparison, we have 
shown that the SXB allows for constraints at lower SFR, 
but that the SXB becomes less sensitive to changes in cx 
at z ^2. Our results in combination with those of Cowie et 
al. (2011) thus provide stronger constraints on the allowed 
redshift evolution of cx. Interestingly, a non-evolution in cx 
with redshift appears at odds with the observed redshift evo- 
lution in the XLFs (unless these are contaminated by low 
luminosity AGN, see above). 

Constraints on the redshift evolution of the Lx-SFR 
relation will be helpful in pinning down the astrophysics that 
is driving the Lx-SFR relation, and may eventually give us 
new insights into the X-ray emissivity of the first galaxies 
which plays a crucial role in determining the the rmal history 
IGM during the dark ages (jMirabel et al.ll201lt ). 
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Figure Al. This plot siiows tiie star formation rate density (co- 
moving) in the Universe, pt(z). The black soli d line shows p* (z) 
that has been derived by iHopkins &: BeacomI ||2006|) . The data 
point at z = represents the more recent z = estimate by Both- 
well et al. (2011). The blue dashed line (red dotted line) shows our 
model in which we evolve i/i* (3>*) in our adopted star formation 
rate function with redshift, while keeping #* (V"*) fixed. Both 
models clearly reproduce the 'observed' redshift evolution of the 
star formation rate density. 
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APPENDIX A: ASSUMED REDSHIFT 
EVOLUTION OF THE STAR FORMATION 
RATE FUNCTION 

Al The Integrated Star Formation Rate Density 

In our paper, we needed to assume the redshift evolution of 
the star formation rate function (dn/dtp). We studied two 
models: 

• In our first model, we evolved ^Ij*{z) to match the ob- 
served star formation rate density, but kept $* fixed. We 
found that the following redshift evolution of ip*{z) provides 
a decent fit to observations: 



riz) 



M, 



yr" 



9.2(1 + 
9.2(1 + 

9.2(1 + zif«fitf 



z < zi; 

Zi ^ z < Z2 

Z > Z2, 



(Al) 



where zi = 1.5 and Z2 = 3.35. The resulting integrated 
star formation rate density is shown as the blue dashed 
line in Figure | A1[ which sh o uld b e compared to that de- 
rived bv I Hopkins fc BeacomI (|2006l ). which is shown as the 
black solid line. Both curves clearly agree. Note that our 
adopted star formation rate function-which was compiled 
from the most recent data-results in a larger star formation 
rate density at 2 = (p, — 0.025A/q yr~^ Mpc"'^, indicated 



by the black filled c ircle a t z = 0), than that inferred by 
iHopkins fc BeacomI (|2006l ). which is pbho6(z = 0) — ah = 
O.OI2M0 yr-i Mpc"^ 

• In our second model, we evolved $*(z) to match the 
observed star formation rate density, but kept tp* fixed. We 
assumed the redshift evolution of ^* (z) was 



$*(2) = $*(2 = 0) X 



PHB06(z) 
PHB06(2 = 0) 



5(2) 



(A2) 



where the function g{z) = ^ + ^(1 -I- z)~^ compensates for 
the fact that the faint-end slope of the star formation rate 
function, a, changes with redshift in our model. The result- 
ing integrated star formation rate density is shown as the 
red dotted line in Figure lAll 



APPENDIX B: DERIVATION OF EQ 1 

Eq [1] plays a central role in our analysis. Here, we pro- 
vide more details on its origin. The total observed flux 
dS from a proper (i.e physical) cosmological volume ele- 
ment dVp is dSiz) = (1 + zfdVpp,£.iz,r)/ATTdl{z). The 
cosmological proper volume element dVp can be written as 



dVp 



Ho {l+z)e(z)"'^P- 



We substitute dAp 



di{z)dn, 



where dAiz) denotes the angular diameter distance to red- 
shift z. We finally get for the differential flux per sterradian 



(Bl) 



dS{z) _ c d\{z)p4z){l + zfC{V,z) _ 
dQ. ^ {1 + z)£{z)-iTvdl{z) 

c pt(z)C(V,z)dz 
i-wHo £{z){l + zy ' 

where we used that dA{z) = (1 + z)~^d-L{z). When we inte- 
grate over redshift and solid angle Ail, we arrive at Eq[T] 



APPENDIX C: CONTRIBUTION OF 
QUIESCENT GALAXIES TO THE SXB 

To compute the contribution of low mass X-ray binaries to 
the SXB we only need to modify Eq[3]in two ways: (i) we re- 
place the star formation rate function with the observed stel- 
lar mass function n(M*,2), and (ii) we replace P{log Lx\tp) 
with P(logLx|M,). The goal of this appendix is to provide 
more details of the calculation, and to show that the conclu- 
sion that LMXBs contribute about an order of magnitude 
less to the SXB than HMXBs is robust to uncertainties in 
the modeling. 

Observed stellar mass functions can be described 
by Schechter functions, and observations have con- 

strai ned the Schechter parame t ers out to z ~ 4 

(e.g. jPerez-Gonzalez et al.f I2OO8I: iMarchesini et all l2009l : 



iKaiisawa et al.ll2009l : lMortlock et al.ll201ll '). In FigurelCTIwe 

show results from calculations in which we adopted the pa- 
rameters from Perez- Gonzalez et al. (2008, red dashed Imes), 
and Mortlock et al. (2011, black solid lines). Both calcula- 
tions agree well. 

We assume that P(logLx|M*) is given by a lognormal 
distribution, wher e (-Lx) = C xA/,. Here, Cx = 8.0 ± 0.5 x 
10^* erg s"^ Mq^ (jGilfanov et al. 2004) . The scatter in this 
relation is not given, and for simplicity we have adopted 
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Figure CI. Same as Figure [21 but for low mass X-ray binaries (LMXB), for which the cumulative luminosity scales linearly with the 
total stellar mass. We integrate the stellar mass functions down to some minimum mass M^^^. The upper left panel shows Sx as a 
function of Mj^in- Another difference with Figure [2] is that we do not extend our calculation beyond Zmax = 4.0 as the stellar mass 
functions are constrained poorly at these redshifts. The black solid lines [red dashed lines] show 5x if we adopt the stellar mass functions 
from Mortlock et al. (2011) [Perez-Gonzalez et al. 2008]. We find that typically, the contribution from LMXBs to the SXB lies about an 
order of magnitude below that of HMXBs. 



(Ti = 0.4, but note that our results can simply be rescaled 
by a factor of exp (^ln'^10[cr2 — ci]) to obtain predictions for 
any (72. The last difference with the calculation presented in 
the main paper is that Lx.lmxb is measur ed in the restframe 
E= 2 - 10 keV band l|Lehmer et al.ll2010l ). 

Figure ICll presents results from our calculations in a 
way that is identical to Figure[2]of the main paper. The main 
differences are: (i) the upper left panel shows Sx as a function 
of minimum stellar mass (instead of minimum star formation 
rate), and (ii) our calculations extend only out to Zmax = 
4.0, as the observed stellar mass functions become uncertain 
there. The general trends in this figure are similar to those 
in Figure[2l except the dependence of Sx on F. This different 
dependence results from the fact that Lx.lmxb is measured 
in the 2-10 keV band (compared to 0.5-8.0 keV for HMXBs) 
which introduces different K-corrections. Generally, we find 
that LMXBs produce Sx ^3 x 10"" erg s"^ cm"^ deg"^, 
which is ~ 10% of the amount contributed by HMXBs. 
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